Search for: All records

Abstract We present a survey of 1D kinetic particle-in-cell simulations of quasi-parallel nonrelativistic shocks to identify the environments favorable for electron acceleration. We explore an unprecedented range of shock speedsv_sh≈ 0.067–0.267c, Alfvén Mach numbers ${\mathcal{M}}_{\mathrm{A}}=5–40$, sonic Mach numbers ${\mathcal{M}}_{\mathrm{s}}=5–160$, as well as the proton-to-electron mass ratiosm_i/m_e= 16–1836. We find that high Alfvén Mach number shocks can channel a large fraction of their kinetic energy into nonthermal particles, self-sustaining magnetic turbulence and acceleration to larger and larger energies. The fraction of injected particles is ≲0.5% for electrons and ≈1% for protons, and the corresponding energy efficiencies are ≲2% and ≈10%, respectively. The extent of the nonthermal tail is sensitive to the Alfvén Mach number; when ${\mathcal{M}}_{\mathrm{A}}\lesssim 10$, the nonthermal electron distribution exhibits minimal growth beyond the average momentum of the downstream thermal protons, independently of the proton-to-electron mass ratio. Acceleration is slow for shocks with low sonic Mach numbers, yet nonthermal electrons still achieve momenta exceeding the downstream thermal proton momentum when the shock Alfvén Mach number is large enough. We provide simulation-based parameterizations of the transition from thermal to nonthermal distribution in the downstream (found at a momentum around $p_{\mathrm{i},\mathrm{e}}/m_{\mathrm{i}}{v}_{\mathrm{sh}}\approx 3\sqrt{m_{\mathrm{i},\mathrm{e}}/m_{\mathrm{i}}}$), as well as the ratio of nonthermal electron to proton number density. The results are applicable to many different environments and are important for modeling shock-powered nonthermal radiation. 

Abstract Pulsars in binary systems with strong companion winds can have the magnetopause separating their magnetosphere from the wind located well within their light cylinder. This bow-like enclosure effectively creates a waveguide that confines the pulsar’s electromagnetic fields and can significantly alter its spindown. In this paper, we study the spindown of compressed pulsar magnetospheres in such systems. We parameterize the confinement as the ratio between the equatorial position of the magnetopause (or standoff distance)R_mand the pulsar’s light cylinderR_LC. Using particle-in-cell simulations, we quantify the pulsar spindown for a range of compressions,R_m/R_LC= 1/3–1, and inclination angles,χ= 0°…90°, between magnetic and rotation axes. Our strongly confined models (R_m/R_LC= 1/3) show two distinct limits. Forχ= 0°, the spindown of a compressed pulsar magnetosphere is enhanced by approximately a factor of three compared to an isolated pulsar due to the increased number of open magnetic field lines. Conversely, forχ= 90°, the compressed pulsar spins down at less than 40% of the rate of an isolated reference pulsar due to the mismatch between the pulsar wind stripe wavelength and the waveguide size. We apply our analysis to the 2.77 s oblique rotator (χ= 60°) in the double-pulsar system PSR J0737-3039. With the numerically derived spindown estimate, we constrain its surface magnetic field toB_*≈ (7.3 ± 0.2) × 10¹¹G. We discuss the time modulation of its period derivative, the effects of compression on its braking index, and implications for the radio eclipse in PSR J0737-3039. 

Abstract Collisionless shocks tend to send charged particles into the upstream, driving electric currents through the plasma. Using kinetic particle-in-cell simulations, we investigate how the background thermal plasma neutralizes such currents in the upstream of quasi-parallel non-relativistic electron–proton shocks. We observe distinct processes in different regions: the far upstream, the shock precursor, and the shock foot. In the far upstream, the current is carried by nonthermal protons, which drive electrostatic modes and produce suprathermal electrons that move toward upstream infinity. Closer to the shock (in the precursor), both the current density and the momentum flux of the beam increase, which leads to electromagnetic streaming instabilities that contribute to the thermalization of suprathermal electrons. At the shock foot, these electrons are exposed to shock-reflected protons, resulting in a two-stream type instability. We analyze these processes and the resulting heating through particle tracking and controlled simulations. In particular, we show that the instability at the shock foot can make the effective thermal speed of electrons comparable to the drift speed of the reflected protons. These findings are important for understanding both the magnetic field amplification and the processes that may lead to the injection of suprathermal electrons into diffusive shock acceleration. 

Abstract Particle-in-cell (PIC) simulations have shown that relativistic collisionless shocks mediated by the Weibel instability accelerate ∼1% of incoming particles, while the majority are transmitted through the shock and become thermalized. The microphysical processes that determine whether an incoming particle will be transmitted or reflected are poorly understood. We study the microphysics of particle reflection in Weibel-mediated shocks by tracking a shell of test particles in a PIC simulation of a shock in pair plasma. We find that electrons in positron-dominated filaments and positrons in electron-dominated filaments efficiently reflect off of strong magnetic structures at the shock. To participate in diffusive shock acceleration, however, these reflected particles headed toward the upstream must avoid getting advected downstream. This is enabled by incoming filaments, which trap reflected particles carrying the same sign of current as the filaments. The final injection efficiency on the order of ∼1% thus results from the effectiveness of the initial reflection at the shock and the reflected particles’ probability of survival in the upstream postreflection. We develop a model that predicts the fraction of high-energy particles as a function of the properties of Weibel filamentation. 

Abstract We study a relativistic collisionless electron–positron shock propagating into an unmagnetized ambient medium using 2D particle-in-cell simulations of unprecedented duration and size. The shock generates intermittent magnetic structures of increasingly larger size as the simulation progresses. Toward the end of our simulation, at around 26,000 plasma times, the magnetic coherence scale approachesλ∼ 100 plasma skin depths, both ahead and behind the shock front. We anticipate a continued growth ofλbeyond the time span of our simulation, as long as the shock accelerates particles to increasingly higher energies. The post-shock field is concentrated in localized patches, which maintain a local magnetic energy fractionε_B∼ 0.1. Particles randomly sampling the downstream fields spend most of their time in low field regions (ε_B≪ 0.1) but emit a large fraction of the synchrotron power in the localized patches with strong fields (ε_B∼ 0.1). Our results have important implications for models of gamma-ray burst afterglows. 

Abstract We present JWST observations of the Crab Nebula, the iconic remnant of the historical SN 1054. The observations include NIRCam and MIRI imaging mosaics plus MIRI/MRS spectra that probe two select locations within the ejecta filaments. We derive a high-resolution map of dust emission and show that the grains are concentrated in the innermost, high-density filaments. These dense filaments coincide with multiple synchrotron bays around the periphery of the Crab's pulsar wind nebula (PWN). We measure synchrotron spectral index changes in small-scale features within the PWN’s torus region, including the well-known knot and wisp structures. The index variations are consistent with Doppler boosting of emission from particles with a broken power-law distribution, providing the first direct evidence that the curvature in the particle injection spectrum is tied to the acceleration mechanism at the termination shock. We detect multiple nickel and iron lines in the ejecta filaments and use photoionization models to derive nickel-to-iron abundance ratios that are a factor of 3–8 higher than the solar ratio. We also find that the previously reported order-of-magnitude higher Ni/Fe values from optical data are consistent with the lower values from JWST when we reanalyze the optical emission using updated atomic data and account for local extinction from dust. We discuss the implications of our results for understanding the nature of the explosion that produced the Crab Nebula and conclude that the observational properties are most consistent with a low-mass Fe core-collapse supernova, even though an electron-capture explosion cannot be ruled out. 

Abstract Some of the most energetic pulsars exhibit rotation-modulatedγ-ray emission in the 0.1–100 GeV band. The luminosity of this emission is typically 0.1%–10% of the pulsar spin-down power (γ-ray efficiency), implying that a significant fraction of the available electromagnetic energy is dissipated in the magnetosphere and reradiated as high-energy photons. To investigate this phenomenon we model a pulsar magnetosphere using 3D particle-in-cell simulations with strong synchrotron cooling. We particularly focus on the dynamics of the equatorial current sheet where magnetic reconnection and energy dissipation take place. Our simulations demonstrate that a fraction of the spin-down power dissipated in the magnetospheric current sheet is controlled by the rate of magnetic reconnection at microphysical plasma scales and only depends on the pulsar inclination angle. We demonstrate that the maximum energy and the distribution function of accelerated pairs is controlled by the available magnetic energy per particle near the current sheet, the magnetization parameter. The shape and the extent of the plasma distribution is imprinted in the observed synchrotron emission, in particular, in the peak and the cutoff of the observed spectrum. We study how the strength of synchrotron cooling affects the observed variety of spectral shapes. Our conclusions naturally explain why pulsars with higher spin-down power have wider spectral shapes and, as a result, lowerγ-ray efficiency. 

Using the field–particle correlation technique, we examine the particle energization in a three-dimensional (one spatial dimension and two velocity dimensions; 1D-2V) continuum Vlasov–Maxwell simulation of a perpendicular magnetized collisionless shock. The combination of the field–particle correlation technique with the high-fidelity representation of the particle distribution function provided by a direct discretization of the Vlasov equation allows us to ascertain the details of the exchange of energy between the electromagnetic fields and the particles in phase space. We identify the velocity-space signatures of shock-drift acceleration of the ions and adiabatic heating of the electrons arising from the perpendicular collisionless shock by constructing a simplified model with the minimum ingredients necessary to produce the observed energization signatures in the self-consistent Vlasov–Maxwell simulation. We are thus able to completely characterize the energy transfer in the perpendicular collisionless shock considered here and provide predictions for the application of the field–particle correlation technique to spacecraft measurements of collisionless shocks.